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ABSTRACT.  A technique  is  proposed  for  the  extraction  of  spectral 
signals  for  which  the  components  within  a signal  set  are  charac- 
terized by  complex  envelopes  which  exhibit  correlation.  The 
technique  is  based  on  an  orthogonal  decomposition  of  the  cross 
(discrete)  frequency  correlation  matrix  describing  the  set  in  the 
presence  of  additive  uncorrelated  noise.  Certain  advantages  of 
the  technique  relative  to  those  based  on  Fourier  analysis  with 
either  estimated  power  spectrum  or  minimum  variance  processing 
are  presented.  Specifically,  it  is  illustrated  how  the  detection 
performance  of  the  technique  is  established  by  the  total  energy 
in  the  spectral  set  rather  than  by  the  levels  of  individual  spec- 
tral elements  in  the  set.  Furthermore,  it  is  shown  how  minimum 
variance  spectrum  analysis  actually  can  suppress  components  in 
such  a set. 

I.  INTRODUCTION 


In  many  spectral  signal  detection  systems  the  detector  imple- 
mentation embodies  a model  for  the  signal  which  assumes  that  the 
signal  is  characterized  by  a random  phase  parameter.  Furthermore, 
it  is  typically  assumed  that  the  behavior  of  amplitude  and  phase 
of  one  spectral  component  is  uncorrelated  with  the  amplitude  and 
phase  of  a spectral  component  which  is  occurring  simultaneously 
in  another  region  of  the  analysis  band.  As  such,  phase  information 
is  ignored  and  detection  is  performed  on  the  basis  of  estimated 
power  spectrum.  More  recently  [1,  2,  3]  sequential  adaptive 
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schemes  have  been  considered  which  exploit  stationary  phase  esti- 
mators which  yield  coherent  detection  algorithms  that  exhibit 
improved  performance.  However,  very  little  consideration  has 
been  given  to  exploiting  potential  correlation  between  the  com- 
plex envelopes  of  multiple  spectral  components  with  the  objective 
of  improving  the  multiple  component  extraction  process.  Intui- 
tively, one  feels  that  the  best  set  extraction  scheme  should 
perform  simultaneous  coherent  detection  and  component  relating. 

In  the  following  section,  such  an  approach  to  coherent  component 
relation  which  occurs  simultaneous  with  detection  is  presented. 
The  proposed  scheme  is  based  on  a principle  component  analysis 
of  the  random  Fourier  transformed  data  vector  (for  example  see 
[4]).  Specifically,  an  orthogonal  decomposition  of  the  cross 
frequency  correlation  (CFC)  matrix  into  its  eigenvectors  indi- 
cates that  the  frequency  information  on  spectral  components  with 
coherent  envelopes  appears  in  a single  eigenvector.  Furthermore, 
the  corresponding  eigenvalue  contains  energy  from  all  spectral 
components.  Accordingly,  it  is  proposed  that  spectral  sets  be 
extracted  by  identification  of  the  appropriate  eigenvalues  and 
eigenvectors. 

II.  THE  CROSS  FREQUENCY  CORRELATION  (CFC)  MATRIX 

Let  the  deterministic  complex  N-dimensional  vector  Ej.  define 
the  position  (frequency)  of  a spectral  component  in  an  N- 
dimensional  complex  sample  space.  A random  premultiplier  s^, 
termed  the  complex  envelope,  accounts  for  identical  multiplica- 
tive random  amplitude  disturbances  and  additive  phase  perturba- 
tions on  each  element  of  the  position  vector  E^.  Consider  an 
observed  N-dimensional  random  data  vector  consisting  of  K < N 
spectral  components  and  additive  noise 

K 

X = £ s.  E.  ♦ N (1) 

k=l 


where  N is  a sample  vector  from  an  identically  distributed  and 
statistically  independent  vector  random  process  with  zero  mean 
and  variance  a^.  If  the  matrix  E with  kth  column  Ej.  and  the 
matrix  P with  l-jth  element  p. . = sTsT  are  introduced,  then  the 
CFC  matrix  ^ 1 ^ 


X X* 


" EPE*  * °0  *N 


(2) 

(3) 


is  defined  wherein  ( )*  is  the  matrix  complex  conjugate  transpose 
operator  and  1^  is  an  N by  N identity  matrix.  Assuming  that  the 
vectors  E^  are  linearly  independent  a matrix  4>  can  be  found  such 
that 
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E = M”1  (4) 

where  4»  has  the  property  <!>**  = I„  and  B is  a K-dimensional  upper 

-1-1* 

triangular  matrix.  The  matrix  B P B is  now  written  in 
diagonal  form  as 

B'1  P B-1*  = Y[y*  (5) 

where  y is  a K by  K matrix  consisting  of  the  orthonormal  eigen- 
-1  -1*  r 

vectors  of  B PB  and  £ is  a K by  K diagonal  matrix  consisting 
of  the  eigenvalues  a2  > a2  > •••  cr2  of  B'1PB-1*. 

Using  (4)  and  (5)  in  (3)  and  post-multiplying  by  4>y  yields 


R4>Y  = 4>Y^k 

(6) 

where 

XK  s £ + °o  h 

(7) 

and 

M - <J>Y 

• 

(8) 
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(9) 

give  the  first  K eigenvalues  and  corresponding  eigenvectors  of 
the  CFC  matrix.  All  remaining  (N-K)  eigenvalues  are  identically 
Oq.  Equation  (9)  indicates  that  the  eigenvectors  of  R (col- 
umns of  M)  are  linear  combinations  of  the  spectral  position  vec- 
tors Ej^  where  the  combination  coefficients  are  functions  of 
(a)  the  inner  products  of  the  position  vectors  E|Ej  through  B 
and  (b)  the  complex  envelope  factor  correlation  characteristics 
through  y. 


To  realize  the  significance  of  (7)  and  (9)  consider  the  case 
of  K spectral  components  with  orthogonal  position  vectors,  i.e. 


Ef  E.  * G.  {..  (6..  = dirac  delta) 

-l  -j  l ij  v ij 

and  complex  envelope  power 

(10) 

Si  * |s.|2. 
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Let  the  K spectral  components  consist  of  L sets  where  the  £th  set 
{i  € contains  components  for  which  the  magnitude  squared  co- 
herence 
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\ 0 for  i € n and  j « jT. 


that  is,  envelope  coherence  exists  only  for  components  within  a 
set.  It  can  be  shown  from  (7)  and  (9)  that 

*A  = E Gi  Si  + 

i€  1 1 0 


E (G.  S.)*  E.  for  £ < L 
i€«A  1 1 _1 


xt  * "o 


for  N-L  < A < N . 


Thus,  detecting  the  Ath  spectral  set  is  accomplished  by  comparing 
an  estimated  to  a2-  and  relating  components  within  the  Ath 
spectral  set  is  accomplished  by  interrogating  to  determine 
which  components  of  the  set  are  present  in  the  linear 

combination  M^.  The  above  case  pertains  to  the  situation  where 
each  of  the  spectral  sets  is  disjoint,  i.e.,  n = 0 with 
A£k.  For  the  case  n / 0,  which  includes  partially  coherent 
complex  envelopes,  leakage  between  sets  occurs.  However,  inclu- 
sion of  a spectral  component  within  a given  set  will  tend  to  be 
in  proportion  to  the  amount  of  coherence  which  exists  between 
that  component  and  the  components  within  the  set. 

III.  THE  RELATION  BETWEEN  ORTHOGONAL  CROSS-FREQUENCY-CORRELATION 
(CFC)  MATRIX  DECOMPOSITION  AND  MINIMUM  VARIANCE  SPECTRUM 
ANALYSIS 

A minimum  variance  estimate,  y^,  of  the  complex  envelope 

for  a possible  component  of  the  random  data  vector  )C  with  spectral 
position  vector  D requires  the  filter 

= (G/D*  R’1  D)  R_1D  (15) 

where  D*I)  = G (see  [7]  for  example) . The  expected  output  power 
for  the  minimum  variance  filter  is 


•yMV'  * (G  R £)  . (] 

whereas  the  expected  output  power  for  a conventional  filtering 
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operation  is 

|yc|2  * D*RD  . (] 

In  the  preceeding  section  it  was  shown  that  the  CFC  matrix  can 
be  written  as 


R ’ £ °\  "A  * «o  !N  • 

Using  (18)  in  (16)  and  (17)  gives 


lyMv!2  * °2  °2/(G  " £ *.l£*M,|2) 

0 1=1  * ~ 


l/cl2  - £ °l  la*«a.l2  * °°l  <2 

2 2 2 

where  = 0$/(°£  + Oq).  Consider  the  special  case  of  two  spec 

tral  components  for  both  coherent  and  incoherent  envelopes. 
Given  orthogonal  position  vectors  with  D*£.  = a.  exp(j<f>.)  such 
that  for  real  cc  (0  < a.  < G)  there  results  from  (19)  1 
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(SNR)V 


(SNR)2a‘ 


1 + (SNR) j 1 + (SNR) 2 


1 + (SNR)X  + (SNR) 2 


[ (SNR) ^a2  ♦ (SNR) 2 a2  + 

2 (SNR)  J (SNR)**  COS(*j-^ 


for  |C, 


and  from  (20) 
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Wr 


= < 


G a2Q  [ (SNR) j aj  + (SNR)2  a \ + 1]  for  |C12|2  = 0 


Ga2  [(SNR^a2  + (SNR)2a2  + 


2 (SNR)  j (SNR)  COSf^-^)  ♦ 1] 


for  iCi2 1 2 = 1 


(22) 


where  (SNR)^  = GS^/o^  is  the  post-filter  signal-to-noise  ratio 
for  £ = with  only  the  ith  spectral  component  present.  It  is 
noted  from  (21)  for  |Cj2|2  = 1 that  the  filter  output  power  for 
estimating  s1  (a^  = G and  a2  = 0)  is  reduced  as  (SNR)2  in- 
creases. This  represents  a suppression  of  the  filtered  signal 
even  though  £ is  well  matched  to  £j.  This  is  because  the  minimum 
variance  filter  assumes  that  s^  ana  s2  are  uncorrelated.  This 
can  be  seen  by  considering  the  sensitivity  expression 
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(SNR) L ♦ 
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(SNR)2][GOq(SNR)1]/ 


aJ)(SNR)1  ♦ 


(SNR)2]‘ 


(23) 


which  is  a monotonical ly  decreasing  function  of  (SNR)2>  No  such 
signal  suppression  effect  is  present  in  either  conventional  fil- 
tering or  signal  set  extraction  using  eigenvalue(-vector)  com- 
ponent relating. 

IV.  CONCLUSION 

A spectrum  analysis  technique  based  on  an  orthogonal  decom- 
position of  the  cross  frequency  correlation  matrix  has  been 
proposed  and  compared  to  both  conventional  linear  and  minimum 
variance  spectrum  analysis.  Orthogonal  component  spectrum  analy- 
sis exhibits  a capability  for  simultaneous  extraction  of  all 
spectrally  disjoint  components  characterized  by  coherence  between 
the  envelopes  of  the  spectrally  disjoint  components.  It  has  been 
shown  that  minimum  variance  spectrum  analysis  for  this  class  of 
signals  can  actually  lead  to  worse  performance  than  conventional 
analysis.  Orthogonal  decomposition,  on  the  other  hand,  exhibits 
detection  performance  governed  by  the  ratio  of  the  total  signal 
power  in  an  envelope  coherent  spectral  set  to  the  noise  level  in 
a single  spectral  resolution  cell.  As  such,  orthogonal  spectrum 
analysis  would  have  application  to  the  extraction  of  a limited. 


I 


I 


but  perhaps  interesting,  class  of  signals  in  both  frequency 
and  spatial  wavenumber  analysis  [5,  6].  Finally,  adaptive 
algorithms  based  on  gradient  search  techniques  are  available 
for  realizing  the  orthogonal  component  spectrum  analysis  tech- 
nique described  herein  [7].  These  algorithms  will  circumvent 
the  necessity  of  first  estimating  and  storing  the  CFC  ma- 
trix with  subsequent  orthogonal  decomposition. 
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